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ABSTRACT 

We present solutions to the accretion disk structure equations in which 
advective cooling is accounted for in a self-consistent way. 

It is shown that for high rates of accretion, above a critical value, when 
without advection there are no solutions extending continuously from large to 
small radii, there exist global solutions with self- consistent advective cooling. 
These solutions are quite different from those solutions obtained earlier by other 
authors, where advective cooling was not treated self-consistently. 

Subject headings: accretion, accretion disks - black hole physics 
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1. Introduction 

Since the end of the seventies it has been realized that for high rates of accretion, M, the 
advection of energy with the accretion flow can play a crucial role in the inner parts of an 
accretion disk around black holes. 

In the last few years a set of papers were published, where it was shown that neglecting 
the advective heat transport at high M leads to qualitatively wrong conclusions about 
the topology of the family of solutions to the system of disk structure equations (see for 
example Abramowicz et al 1995; Chen et al 1995; Artemova et al 1996). 

The disk structure equations without advection give rise two branches of solutions: 
optically thick and optically thin, which do not intersect for L < L b < (0.6 — 0.9)L Edd for 
a = 1 and M BH = 10 8 M Q , where L EAd is the Eddington luminosity (Artemova et al 1996). 
For larger luminosities there are no solutions of these equations extending continuously 
from large to small radii. 

It was argued by Artemova et al (1996), that for luminosities larger than Lf,, advection 
becomes important and induces a transition from the optically thick to the optically thin 
disk in the global solution. Thus, if we were to take advection into account, solutions to the 
disk structure equations should also exist at L > L b . In that case the disk should consist of 
an optically thin inner part, smoothly connected to an optically thick outer part. 

The goal of this letter is to construct explicitly accretion disk models taking advective 
heat transport into account self- consistently and to demonstrate in such a way the 
correctness of our hypothesis, as proposed in the paper by Artemova et al (1996). 

We use from now on geometric units with G — 1, c = 1, and use r as the radial 
coordinate scaled to r g = GM/c 2 . 
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2. The Model and the Method of Solution 

In describing accretion disks when advective cooling is important, one can use the 
following model: We assume that advective cooling is added in a self- consistent way to 
other (local) cooling mechanisms in a Keplerian, geometrically thin disk. We neglect the 
radial momentum equation and suggest that the disk thickness is everywhere small. This 
simplification is not too serious because neglecting these factors leads to quantitative 
corrections only, while neglecting advection cooling leads to qualitatively different solutions. 
This was discussed in detail by Abramowicz et al (1995) and Chen et al (1995). Their 
treatment of advective cooling was not self- consistent, however, as their advection 
parameter, £, was assumed to be constant throughout the disk as well as being independent 
of the accretion rate. 

Thus we assume that the disk is geometrically thin and the rotational velocity 
is Keplerian under the pseudo-Newtonian potential proposed by Paczyhski and Wiita 
(1980), $ = —M/(r — 2). We normalize the accretion rate asm = M/M Edd , where 
M Edd = L Edd /c 2 = 4nMm p /a T . 

The non-relativistic disk structure equations without advection, in the case of the 
Newtonian potential, together with the method of their solutions have been described 
earlier (see Artemova et al 1996). Here we will outline the main difference related to the 
use of the potential of Paczyhski and Wiita and to the inclusion of the advective cooling 
process in the energy balance. Paczyhski and Bisnovatyi-Kogan (1981) were the first to 
describe a method to account for the advective cooling in the disk accretion problem. 

The angular velocity for the Paczyhski and Wiita potential is il = [M/r(r — 2) 2 ] 1 / 2 . 
Recall that conservation of angular momentum for a steady-state accretion in the a-disk 
model can be written as 
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where the factor, / = (1 — r? n fi(r in )/(r 2 f2(r))), accounts for the boundary condition at the 
inner edge of the disk, r in . Neglecting advection we have r in = 6. The half thickness of the 
disk is denoted by h(r), and P(r) is the total pressure in the equatorial plane of the disk. 

The viscous heating rate per unit area, Q + , is given by the formula (see e.g. 
Bisnovatyi-Kogan, 1989; Frank, King & Raine, 1992) 
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The advection cooling rate can be written in the form (see Chen and Taam 1993): 
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where T(r) is the temperature, S(r) is the specific entropy, and ^(r) is a factor which 
characterizes the advective cooling. This factor can be written in the form: 
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where, E is the energy per unit mass of the gas, p is the matter density and v = 1/p. The 
energy equation has the form Q + = Q a dv + Q\oc, where Q\ oc is the rate of all local cooling 
processes (see Artemova et al 1996). 

All other equations are the same as in Artemova et al (1996). The principal difference 
here is that now the energy balance is described by a differential equation rather than an 
algebraic one. 

The boundary conditions for this system are the following: at large, r » 100r 9 , the 
solution must coincide with the solution obtained without advection, and at the inner disk 
boundary we should not require that r in = 6, but rather we treat r in as an eigenvalue of the 
problem and allow it to be less than 6. 
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We solved this system of equations by the method of subsequent iterations, with fixed 
rh and a, and varying r in to obtain a self- consistent solution. We did this in two steps: 
started one calculation from large radii and another at the inner edge of the disk and 
adjusting r in until the solutions matched near r R3 13. 

3. Results and Discussion 

We will compare the solutions with and without advection cooling. 

In Figure 1 we plot the optical depth r (oc £, the disk surface density) as a function 
of radius, r, for the case M BH = 1O 8 M , a — 1; for various values of rh, without taking 
into account the advection cooling, i.e. we set £ = 0. For accretion rates rh < fn b = 14.315, 
(corresponding to L 0.9L Edd ), there are two families of solutions which are continuous 
from large radii to the inner disk edge. For rh > 14.315 there are no solutions for a range 
of radii around r R3 13, thus for such high accretion rates there are no solutions extending 
continuously from large to small radii. (See also the discussion in the paper Artemova et 
al, 1966). 

Let us now include advection cooling. For rh < 13 this gives rather small corrections 
to the solutions without advection. When rh > 13 advection is essential. 

In Figure 2 we compare the solutions for rh = 14.0 and rh = 15.0 without advection, 
(£ = 0), and with self-consistent advection, £ = £(r, rh). The functions £(r, rh) for the cases 
rh\ = 14.0 and rh 2 = 15.0 are shown in Figure 3. The eigenvalues of r in are correspondingly: 

r in,riu ~ 6.0 J r in m 2 ~ 5.4. 

In Figure 4 we plot the distribution of pressure P = P(r) for rh\ = 14.0 and m 2 = 15.0 
In the case rh\ = 14.0 the solution changes quantitatively, but not qualitatively. The 
solution is optically thick everywhere in the disk. 
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For the value m = 15.0 > rh b = 14.315 there was no solution extending continuously 
from large to small radii without advection, i.e. when £ = (see Fig.l). With the 
self- consistent advection term £ = £(r) (see Fig. 3) there is a solution which consists of an 
optically thin inner part that continuously turns into an optically thick outer part (see Fig. 

2) . Self- consistent solutions in this case of course also exist for higher values of m. 

Thus, the results presented here demonstrate that the hypothesis proposed in the 
paper by Artemova et al (1996) is correct, namely that for m > m b there are solutions to 
the disk structure equations with self-consistent advective cooling and with transitions from 
optically thick to optically thin regions. Schemes showing the topology of the families of 
disk solutions without advection or with self-inconsistent advection (£ = const, for example) 
should be revised. The advection term, £(r, m), crucially depends on r and m (recall Fig. 

3) . The opposite conclusion that £ does not depend very much on rh was obtained in the 
paper by Chen (1995). 

We note also that there are solutions for m > m b with self- consistent £(r), with 
optically thick inner part and an optically thin outer region of the disk. We will address 
these solutions as well as the more general solution, including the radial pressure gradient 
and angular velocities of the gas, in future work. 

This paper was supported in part by the Danish Natural Science Research Council 
through grant 11-9640-1, in part by Danmarks Grundforskningsfond through its support for 
the establishment of the Theoretical Astrophysical Center. G. B. thanks NORDITA and 
TAC, and G.B.-K. thanks TAC for partial support and hospitality during the initial stages 
of this work. 
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FIGURE CAPTIONS 

Fig. 1. — Optical depth tq, (oc So) as a function of radius, r, for Mbh = 1O 8 M , a = 1.0, 
without advective cooling (£ = 0) and different values of rh. The solid curves correspond 
to optically thin (lower branch) and optically thick (upper) solutions for rh = 14.0. The 
dashed curves are calculated for rh = 15.0 > rh b . No solution extending continuously from 
the outer disk to its innermost edge exists in the latter case. 

Fig. 2. — Comparing the radial dependence of the optical depth r , (oc £ ) for models 
with the same parameters as in Fig. 1. We show the solutions without advective cooling 
(£ = 0) and with self-consistent advective cooling. The corresponding £ = £(r, rh) are shown 
in Fig. 3. The solid curves corresponds to rh = 14.0, £ = 0, the dotted curve corresponds to 
rh = 14.0, £ = £(r, rh), the dashed curves correspond to rh = 15.0, £ = 0. The dot-dashed 
and dashed-triple dotted curve correspond to rh — 15.0, £ = £(r, rh = 15.0). 

Fig. 3. — The quantity £ as a function of radius, r, and the accretion rate rh, for 
the same disk parameters as in Fig. 1. The dotted curve corresponds to rh — 14.0, the 
dot-dashed and dashed-triple dotted curves correspond to rh = 15.0. 

Fig. 4. — The pressure structure of the same disk models as in Fig. 2, (in units 
V = P/pc 2 ). The styles of the curves are the same as in Fig. 2. 
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